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Abstract 

Since the quark-gluon plasma (QGP) reveals some obvious similarities to the well-known electro- 
magnetic plasma (EMP), an accumulated knowledge on EMP can be used in the QGP studies. 
After discussing similarities and differences of the two systems, we present theoretical tools which 
are used to describe the plasmas. The tools include: kinetic theory, hydrodynamic approach 
and diagrammatic perturbative methods. We consider collective phenomena in the plasma with 
a particular emphasis on instabilities which crucially influence temporal evolution of the system. 
Finally, properties of strongly coupled plasma are discussed. 
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1 INTRODUCTION 



The plasma - the ionized gas of electrons and ions - has been actively studied 
since its discovery in a discharge tube at the end of 19-th century. The term 
plasma was introduced by Irving Langmuir in 1929. Prospects to get a practically 
unlimited source of energy due to nuclear fusion reactions in a hot ionized gas 
of hydrogen isotopes, have stimulated a large scale program to study plasmas in 
terrestrial experiments for over half a century. Plasmas are also actively studied 
by astrophysicists as it appears the most common phase of matter. About 99% 
of the entire visible Universe is in the plasma phase, not only stars are formed 
of ionized gas but the interstellar and intergalactic medium is also a plasma, 
although a very sparse one. Principles of the plasma physics can be found in e.g. 
the well-known textbooks {H[2]). 

The quark-gluon plasma (QGP) is the system of quarks and gluons which 
are not confined in hadron's interiors but can freely move in a whole volume 
occupied by the system. A broad presentation of the whole field of QGP physics 
is contained in three volumes of review articles (OHHU); the lectures ([6]) can 
serve as an elementary introduction. Active studies of QGP started in the mid 
1980s when relativistic heavy-ion collisions offered an opportunity to create a 
drop of QGP in a laboratory. The experimental programs at CERN and BNL 
provided an evidence of the QGP production at the early stage of nucleus-nucleus 
collisions, when the system is extremely hot and dense, but properties of QGP 
remain enigmatic. So, one can ask: what do electromagnetic plasmas (EMP) 
tells us about QGP? 

The QGP reveals some obvious similarities to the well-known electromagnetic 
plasma, as Quantum Chromodynamics (QCD) describing the interactions of the 
quarks and gluons resembles Quantum Electrodynamics (QED) which governs 
interactions of charged objects. Thus, some lessons from EMP should be useful in 
the exploration of QGP. The aim of this article is to discuss what QGP physicists 
can actually learn from their EMP colleagues and how the huge accumulated 
knowledge on EMP can be used in the QGP studies. However, we have to be 
aware not only of similarities but of important differences between EMP and 
QGP. Some differences are of rather trivial origin but some are deeply rooted in 
dynamical foundations of the two systems. 

Let us enumerate these 'trivial' dissimilarities. QGP is usually relativistic or 
even ultrarelativistic while EMP is mostly nonrelativistic in laboratory experi- 
ments. The differences between the non-relativistic and relativistic plasmas go far 
beyond the kinematics of motion of plasma particles. For example, let us consider 
the plasma's composition. In the non-relativistic system, there are particles but 
no antiparticles and the particle's number is conserved. In the relativistic system 
we have both particles andantiparticles, as electrons and positrons in EMP, and 
the lepton number - not the particle's number - is conserved. Particle number 
density is not a right quantity to characterize the system. For this reason, the 
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QGP physicists use the baryon and strangeness densities. 

Another 'trivial' but very important distinctive feature of EMP is the huge 
mass difference between electrons and ions which is responsible for a specific dy- 
namic role of heavy ions. The ions are usually treated as a passive background, 
which merely compensates the charge of electrons, but electro- ion collisions drive 
the system towards equilibrium and maintain the equilibrium. However, the 
energy transfer between electrons and ions is very inefficient and their mutual 
equilibration is very slow. Therefore, we have electron and ion fluids of different 
temperatures for a relatively long time. There is nothing similar in QGP. There 
are heavy quarks - charm, bottom and top - which are, however, much less popu- 
lated than the light quarks and gluons and their lifetime is short. Therefore, the 
heavy quarks hardly influence the QGP dynamics. 

The electromagnetic plasma, which is the closest analogue of QGP, is the rel- 
ativistic system of electrons, positrons and photons. Such a plasma is actually 
studied in context of some astrophysical applications, e.g. supernovae explosions. 
Then, the differences between QGP and EMP are of dynamical origin - the first 
one is governed by Quantum Chromodynamics while the second one by Quantum 
Electrodynamics. The former theory is Abelian while the later one is nonAbelian 
with a prominent role of gluons which carry color charges and thus, not only me- 
diate the interaction among colored quarks and antiquarks but interact among 
themselves. Gluons, in contrast to photons, also contribute to the density of color 
charges and to the color current. 

The most important common feature of EMP and QGP is the collective char- 
acter of the dynamics. The range of electrostatic interaction is, in spite of the 
screening, usually much larger than the inter-particle spacing. There are many 
particles in the Debye sphere - the sphere of the radius equals the effective in- 
teraction range, and motion of these particles is highly correlated. There is a 
similar situation in the deconfined perturbative phase of QCD ([7]). The Debye 
mass is of order gT where g is the QCD constant and T is the temperature. Since 
the particle density in QGP is of order T 3 , the number of partons in the Debye 
sphere, which is roughly 1/g 3 , is large in the weakly coupled (1/g S> 1) QGP. 

In various laboratory experiments, EMP is embedded in an external electro- 
magnetic field. For example, the magnetic field is used to trap the plasma, and 
there are numerous fascinating phenomena occurring in such a situation. In the 
case of QGP produced in relativistic heavy-ion collisions, it is hard to imagine any 
external chromodynamic field applied to the plasma. Therefore, we will consider 
here only the systems where fields are self-consistently generated in the plasma. 

Our article is organized as follows. Theoretical tools, which are used to describe 
the plasmas, are presented in Chapter[2j The tools include: the kinetic theory, the 
hydrodynamic approach and diagrammatic methods of field theory. In Chapter [3] 
we discuss collective phenomena which are the most characteristic feature of 
plasmas. After explaining the phenomenon of screening, quasi-particle modes 
in the equilibrium and non-equilibrium plasma are presented. We pay much 
attention to instabilities which crucially influence plasma dynamics. The problem 
of particle's energy loss in a plasma is also discussed. Chapter |4] is devoted to 
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the strongly coupled plasma which reveals particularly interesting properties. 

Throughout the paper we use the natural units with c = H = ks = 1 and 
the metric (1, — 1, — 1, — 1). However, there is a little complication here. Plasma 
physicists usually use the Gauss (CGS) units, where the fine structure constant 
equals a = e 2 ~ 1/137, while the electromagnetic counterpart of the units usually 
applied in QCD is the so-called Heaviside-Lorentz system where the 4ir factor 
does not show up in the Maxwell equations but a = e 2 /4-7r. We stick to the 
traditionally used units in the two fields of physics, and thus the factor of 4-7T has 
to be additionally taken into account to compare EMP to QGP formulas. 

2 THEORETICAL TOOLS 
2.1 Transport theory 

Transport theory provides a natural framework to study equilibrium and nonequi- 
librium plasmas. The central object of the theory is the distribution func- 
tion which describes a time dependent distribution of particles in a phase-space 
spanned by the particle's momenta and positions. The distribution function of 
each plasma component evolves due to the inter-particle collisions and the inter- 
action with an external and/or self-consistently generated mean-field. The two 
dynamical effects give rise to the collision and mean-field terms of a transport 
equation satisfied by the distribution function. 

2.1.1 ELECTROMAGNETIC PLASMA A formulation of the kinetic 
theory of relativistic plasma can be found in (jSJ). The distribution function is 
denoted as f n (p,x) with the index n labeling plasma components: electrons, 
positrons, ions. Spin is usually treated as an internal degree of freedom. The 
function depends on the four-position x = (t, x) and the three- momentum p. 
The four- momentum p obeys the mass-shell constraint p 2 = m 2 , where m is the 
particle mass, and then p = (E p , p) with E p = \J m 2 + p 2 . 
The distribution function satisfies the transport equation 

(p^ + QnP^F^d;) f n (p, X) = C[f n ) , (1) 

where C[f n ] denotes the collision term, q n is the charge of the plasma species 
n and F^ u is the electromagnetic strength tensor which either represents an 
external field applied to the system or/and is generated self-consistently by the 
four-currents present in the plasma 

d^F^ = 4tt> , 

where 

The transport equation can be solved in the linear response approximation. 
The equation is linearized around the stationary and homogeneous state described 
by the distribution f n (p)- The state is also assumed to be neutral and there are 
no currents. The distribution function is then decomposed as 

fn(p, x) = fn(p) + Sf n (p, x) , 
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where / n (p) > 6f n (p,x). 

The transport equation linearized in 5f n and F^ u can be exactly solved after 
the Fourier transformation, which is defined as 

/(*) = J otx e ik *f(x) , f(x) = / ^4 e~ ikx f(k) . (3) 

Then, one finds 5f n (p,k), which is the Fourier transform of 5f n (p,x), and the 
induced current, which can be written as 

Sf(k)=-U^(k)A l/ (k), (4) 
with the polarization tensor equal to 

d 3 p - (p • k) 2 g» v + fc 2 pV _ (p . fc)(jfe/y + fc" p M) 



{2nY Jn ^> (p-k) 2 

(5) 

The tensor is symmetric (H^ u (k) = IL VfJ, (k)) and transverse (k^H^ u (k) = 0) which 
guarantees that the current given by Equation [J] is gauge independent. 

For isotropic plasmas, the polarization tensor has only two independent com- 
ponents which are usually chosen as 

U L (k) = Iloo(fc) , 

Ht(*) = \ (^-^) (6) 

where the indices = 1,2,3 label three-vector and tensor components. In the 
case of an ultrarelativistic (T 3> to) electron-positron equilibrium plasma, the 
momentum integral in Equation [5] can be performed analytically in the high- 
temperature limit (T 3> lu, |k|), and the result derived already in 1960 by Silin 
(El) reads 




n L (fc) 



(7) 



where A; = k) and to 7 = eT/3 denotes the thermal photon mass generated by 
the interaction of the photons with the electrons and positrons. 

The above polarization tensor was found in the collisionless limit of the trans- 
port equation. The effect of collisions can be easily taken into account if the 
so-called BGK collision term is used in the transport equation (flQ|) . The result 
for an ultrarelativistic equilibrium plasma is given in pip . 

2.1.2 QUARK-GLUON PLASMA The transport theory of QGP (pjfTBl) 
appears to be much more complicated than its electromagnetic counterpart. The 
distribution function of quarks Q(p, x) is a hermitian N c xN c matrix in color space 
(for a SU(iV c ) color gauge group). The distribution function is gauge dependent 
and it transforms under a local gauge transformation U(x) as 

Q(p,x) -+U(x)Q(p,x)rt(x) . (8) 
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Here and in the most cases below, the color indices are suppressed. The distri- 
bution function of antiquarks, which we denote by Q(p,x), is also a hermitian 
N c x N c matrix and it transforms according to Equation [8l The distribution 
function of gluons is a hermitian (JV| — 1) x — 1) matrix which transforms as 

G(p,x) -^U{x) G{p,x)U\x) , (9) 

where 

U ab (x) =2Ti[T a U{x)r h U\x)] , 

with r a , a = 1, ...,N^ — 1 being the SU(iV c ) group generators in the fundamental 
representation with Tr(r a r b ) = \5 ab . 

The color current is expressed in the fundamental representation as 

f(?) = - f / (03 ^ [QM-QM (io) 

- ^-Tr[Q(p, x) - Q(p, x)] + 2r a Tr[TC(p, x)] 

where g is the QCD coupling constant. A sum over helicities, two per particle, 
and over quark flavors Nf is understood in Equation [TUl even though it is not 
explicitly written down. The SXJ(N C ) generators in the adjoint representation 
are expressed through the structure constants T^ c = —if a bc, and are normalized 
as Tr [T a T b ] = N c 5 ab . The current can be decomposed as j fl (x) = j£(x)T a with 
3a( x ) = 2Tr(r a j' M (x)). The distribution functions, which are proportional to the 
unit matrix in color space, are gauge independent and they provide the color 
current (Equation I10p which identically vanishes. 

Gauge invariant quantities are given by the traces of the distribution functions. 
Thus, the baryon current and the energy-momentum tensor read 



(2tt) e 

T^(x) = |||//Tr[Q(p,x) + Q(p,x) + G(p,x)] ) 

where we use the same symbol Tr[- • •] for the trace both in the fundamental and 
adjoint representations. 

The distribution functions of quarks, antiquarks and gluons satisfy the trans- 
port equations: 

p^D fl Q(p,x) + ^{F tll/ (x),d;Q(p,x)} = C[Q,Q,G], (11) 
^ J D / ,Q(p,x)-|^{^(x),^Q(p,x)} = C[Q,Q,G], (12) 

p^G(p,x) + |^{^(x),^G(p,x)} = C g [Q,Q,G], (13) 

where {..., ...} denotes the anticommutator and dp the four-momentum deriva- 
the covariant derivatives and T>^ act as 

0^ = 8^- ig[A^(x), ... ] , V^ = d^- ig[A^{x), ... ] , 



1 As the distribution functions do not depend on po, the derivative over po is identically zero. 
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with and An being four-potentials in the fundamental and adjoint represen- 
tations, respectively: 



The strength tensor in the fundamental representation is = d^A v — d v A^ — 
ig[A^,A u ], while denotes the field strength tensor in the adjoint representa- 
tion. C, C and C g represent the collision terms. 

The transport equations are supplemented by the Yang- Mills equation describ- 
ing generation of the gauge field 



where the color current is given by Equation [TUJ 

As in the case of the EM plasma, the transport equations, which are linearized 
around a stationary, homogeneous and colorless state can be solved. Because 
of the color neutrality assumption, the analysis is rather similar to that of the 
Abelian plasma, and it ends up with the polarization tensor which is proportional 
to the unit matrix in the color space and has the form of Equation [5j 

As in the case of EMP, the collisions can be easily taken into account using the 
approximate BGK collision terms (|14|ll5p . Within a more realistic approach color 
charges are treated in a similar way as spin degrees of freedom, and ones uses the 
so-called Waldmann-Snider collision terms (fT6|[T7]) which are usually applied to 
study spin transport. 

2.2 Hydrodynamic approach 

Within the hydrodynamic approach, the plasma is treated as a liquid and it is 
described in terms of macroscopic variables which obey the equations of motion 
resulting from the conservation laws. The fluid equations are applied to a large 
variety of plasma phenomena but, depending of the time scale of interest, the 
actual physical content of the equations is rather different. 

Real hydrodynamics deals with systems which are in local equilibrium, and 
thus it is only applicable at sufficiently long time scales. The continuity and the 
Euler or Navier-Stokes equations are supplemented by the equation of state to 
form a complete set of equations. The equations can be derived from kinetic 
theory, using the distribution function of local equilibrium which by definition 
maximizes the entropy density, and thus, it cancels the collision terms of the 
transport equations. 

In the electron-ion plasma there are several time scales of equilibration. The 
electron component of the plasma reaches the equilibrium in the shortest time, 
then ions are equilibrated but for a relatively long time the electron and ion 
temperatures remain different from each other, as the energy transfer between 
electrons to ions is rather inefficient. This happens due to the huge mass difference 
between electrons and ions. 

When the electrons have reached local equilibrium with their own temper- 
ature and hydrodynamic velocity, the collision terms of the kinetic equations 



A»(x) = A»{x)t' 



A»{x) = T a A%(x) . 




(14) 
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representing electron-electron collisions vanish while the collision terms due to 
electron-ion collisions can be neglected as they influence the electron distribution 
function only at a sufficiently long time scale. Then, one obtains hydrodynamic 
equations of an electron fluid. When the ion component is also equilibrated we 
have two fluids with different temperatures and hydrodynamic velocities. At the 
times scales when the fluid equations are applicable, the plasma can be treated 
as locally neutral. Charge fluctuations are obviously possible but they disappear 
fast as the electric field generated by the local charges induces the currents which 
in turn neutralize the charges. Since the plasma is nearly an ideal conductor, the 
process of plasma neutralization is very fast. Due to the charge neutrality of the 
plasma, the electric field is not present in the fluid equations and we end up with 
the magnetohydrodynamics where the pressure gradients and magnetic field drive 
the plasma dynamics. 

As explained above, the regime of magnetohydrodynamics appears because 
there is a heavy positive component of the plasma (ions) and a light negative 
component (electrons). There is no QCD analogue of magnetohydrodynamics 
as every quark or gluon can carry opposite color charges. Therefore, when local 
equilibrium is reached various color components of the plasma have the same 
temperatures and hydrodynamic velocities (fTT|) . Since the quark-gluon system 
becomes color neutral even before the local equilibration is reached (I16pi4p . we 
deal with hydrodynamics of neutral fluid where the chromodynamic fields are 
absent. Such a relativistic hydrodynamics of colorless QGP has been actively 
studied over the last two decades (|18pi9p . 

The hydrodynamic equations, which actually express macroscopic conservation 
laws, hold not only for systems in local equilibrium but for systems out of equi- 
librium as well. The equations can be then applied at time scales significantly 
shorter than that of local equilibration. At such a short time scale, the colli- 
sion terms of the transport equations can be neglected entirely. However, extra 
assumptions are then needed to close the set of equations, as the (equilibrium) 
equation of state cannot be used. Plasma physicists developed several methods 
to close the set of equations, and thus fluid equations are used to study bulk 
features of short time scale phenomena in the plasmas. To get more detailed in- 
formation, the kinetic theory is needed. Since the fluid equations are noticeably 
simpler than the kinetic ones, the hydrodynamic approach is frequently used in 
numerical simulations of plasma evolution, studies of nonlinear dynamics, etc. 

Below we derive the fluid equations for EMP and QGP from the respective 
kinetic theory. Since the fluid approach under consideration is supposed to hold 
at sufficiently short time scales, we use the collisionless transport equations. 
2.2.1 ELECTROMAGNETIC PLASMA We assume here that there are 
several streams in the relativistic plasma system and that the distribution func- 
tions of each plasma component (electrons, positrons, ions) belonging of each 
stream satisfy the collisionless transport equation. The equations are coupled 
only through the electromagnetic mean field which is generated by the current 
coming from all streams. The field in turn interacts with every stream. 

Integrating the collisionless transport Equation Q] over momentum, one finds 
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0X = , (15) 

where the four-flow is 

<{X) = J j^LpPf^x). (16) 

The index a labels simultaneously the streams and plasma components. 

Multiplying the transport Equation Q] by the four-momentum p and integrating 
it over momentum, we get 

^rr + <fa<V = ' ( 17 ) 

where the energy-momentum tensor is 

T£T{x) = J -0^p^f a (p,x). (18) 

The structure of n£ and T£ v is assumed to be that of the ideal fluid in local 
thermodynamic equilibrium. Thus, one has 

<(x) = n Q (x)<(x), (19) 

TH u {x) = [e a (x)+p a (x)]u^x)u v (x)-p a (x)g^ . (20) 

To obtain the relativistic version of the Euler equation, Equation [T7] needs to 
be manipulated following (20). Substituting the energy-momentum tensor of the 
form of Equation [20] into Equation [17] and projecting the result on the direction 
of u£, one finds 

u av d^ u = <d^e Q + (e a + p a )^< = . (21) 

Computing d^T^ —u^Uapd^T^P , one gets the Lorentz covariant form of the Euler 
equation 

M v a = (e a + p a )u afl d"u v a + («d M - 8 v )p - q a n a u afl F^ = . (22) 
The equation in a more familiar form is given by M a - v a M° = 0. Namely, 

(ea+p Q )7a(^+ v Q- v ) v Q + (v+v Q — )p tt -g a Ji a 7 [E-v a (v a 'E)+v a xB] =0 , 

(23) 

where the four-velocity was expressed as = (7 a ,7a v a ) with j a = (1 — 

In the nonrelativistic limit (which is easily obtained when the velocity of light 
c is restored in the equation) , Equation [23] gets the well-known form 

(7* + v « • V ) v « + " — ( E + v Q x B) = . (24) 

Var / m a n a m a 

The fluid Equations HU H7] with and T£ v given by Equations HU [20] do 
not constitute a closed set of equations - there are 5 equations and 6 unknown 
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functions: n a , p a , e a and 3 components of u£ (because of the constraint u^u^a = 
1, one component of u 1 ^ can be eliminated). There are several methods to close the 
set of equations. In particular, assuming that the system's dynamics is dominated 
by the mean- field interaction, one can neglect the pressure gradients. One can 
also add an equation which relates p a to e a . The relation is usually called the 
equation of state, but one should be aware that the plasma system is not in 
equilibrium, and the thermodynamic relations in general do not hold. 

In the ultrarelativistic limit when the characteristic particle's energy (the tem- 
perature of the equilibrium system) is much larger than the particle's mass, and 
thus p 2 = 0, the energy- momentum tensor is traceless (Tjf a = 0), as follows 
from Equation [18] for p 2 = 0. Then, Equation [20] combined with the constraint 
u^(x)u a ^(x) = 1 provides the desired relation 

e a (x) = 3p a (x) , (25) 

which coincides with the equation of state of an ideal gas of massless particles. 

Since the distribution functions of every plasma component belonging to every 
stream are assumed to obey the collisionless transport equation, we have a sep- 
arated set of fluid equations for every plasma component of every stream. The 
equations are coupled only through the electromagnetic mean field. More pre- 
cisely, the electrons, positrons and ions of every stream contribute to the current 
generating the field which in turn interacts with the streams. 

The fluid equations can be solved in the linear response approximation. The 
equations are linearized around the stationary and homogeneous state described 
by n a and n^. The state is neutral and there are no currents i.e. 



E 



n a < = . (26) 



The charge density is decomposed as 

n a (x) = n a + 5n a (x) , (27) 

where n a ^> 5n a . The fully analogous decomposition of the hydrodynamic ve- 
locity pressure p a and energy density e Q is also adopted. 

The set of the continuity and Euler equations linearized in 5n a , du 1 ^, 5p a , 5e a , 
and F^ u can be exactly solved after they are Fourier transformed. Thus, one 
finds 5n a (k), 5u^(k) when the set of fluid equations is closed by neglecting the 
pressure gradients. If the equation of state is used, one also finds 5e a (k). 

Keeping in mind that the induced current equals 

5f = ^2(q a n a 5u£ + q a 5n a u%) , 

a 

one finds from Equation [J] 

W u {k) = jj^jj 1 \k 2 KK + (u Q • fc) V ~ {u a ■ k)(k^ a + feX) 

a + Pa \U a ' ft J L 

(u a ■ k)k 2 (k^u v a + k v u£) - (u a ■ k) 2 k»k v - k A u^u u a 
k 2 + 2(u a ■ k) 2 



+ 



(28) 
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The first term gives the polarization tensor when the pressure gradients are 
neglected while the second term gives the effect of the pressure gradients due 
to the equation of state given by Equation [25] The first term is symmetric 
(Tl^(k) = H l '^{k)) and transverse {k^Yl^^k) = 0). The second term is symmet- 
ric and transverse as well. Thus, the whole polarization tensor (Equation I28[) is 
symmetric and transverse. The first term of Equation [28] can be obtained from 
the kinetic theory result (Equation [5]) with the distribution function / n (p) pro- 
portional to <5( 3 )(p — (e a + p a )u a /n a ). Thus, the first term neglects the thermal 
motion of plasma particles while the second term takes the effect into account. 

2.2.2 QUARK-GLUON PLASMA The fluid approach presented here fol- 
lows the formulation given in (21). As in the EMP case, we assume that there 
are several streams in the plasma system and that the distribution functions of 
quarks, antiquarks and gluons of each stream satisfy the collisionless transport 
equation. The streams are labeled with the index a. 

Further analysis is limited to quarks but inclusion of anti-quarks and gluons 
is straightforward. The distribution function of quarks belonging to the stream 
a is denoted as Q a (p, %)■ Integrating over momentum the collisionless transport 
Equation 1111 satisfied by Q a , one finds the covariant continuity equation 

£>X = » ( 29 ) 
where n£J is N c x iV c matrix defined as 

<(x) = J j^lfQafax). (30) 

The four-flow transforms under gauge transformations as the quark distribu- 
tion function, i.e. according to Equation [8] 

Multiplying the transport Equation Qj] by the four- momentum and integrating 
the product over momentum, we get 

^r-f{v><} = °' ( 31 ) 

where the energy-momentum tensor is 

TST{x) = J -0^p^Q a (p,x). (32) 

We further assume that the structure of and Tg v is 

<(x) = n a (x)u%(x), (33) 

TTix) = ±(e a (x)+p a (x)){uZ(x),u» a (x)}-p a (x)g^ , (34) 

where the hydrodynamic velocity is, as n a , e a and p a , a N c x N c matrix. 
The anticommutator of and u v a is present in Equation [M] to guarantee the 
symmetry of T£ u with respect to fi «-> v which is evident in Equation [321 

In the case of an Abelian plasma, the relativistic version of the Euler equation 
is obtained from Equation [31] by removing from it the part which is parallel to 
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u£. An analogous procedure is not possible for the non-Abelian plasma because 
the matrices n a , and it„, in general, do not commute with each other. Thus, 
one has to work directly with Equations [29l [3T1 with and Tg u defined by 
Equations E31 EH The equations have to be supplemented by the Yang-Mills 
Equation Q3] with the color current of the form 



where only the quark contribution is taken into account. 

The fluid Equations EH EH as their EM counter part, do not form a closed 
set of equations but can be closed analogously. The only difference is that the 
equation of state (Equation |25|) relates to each other the matrix value functions 
e a and p a . 

As in the case of the EM plasma, the fluid Equations ll5l [T71 which are linearized 
around a stationary, homogeneous and colorless state described by h a , e a , p a and 
u^, can be solved (pTj) . Because of the color neutrality assumption, n a , e a , p a and 
are all proportional to the unit matrix in the color space, the analysis is rather 
similar to that of the Abelian plasma, and one ends up with the polarization 
tensor from Equation [28] which is proportional to the unit matrix in the color 
space. 

2.3 Diagrammatic methods 

Various characteristics of the weakly coupled plasma can be calculated using 
the perturbative expansion, that is diagrammatic methods of field theory. It 
requires a generalization of the Feynman rules applicable to processes, which 
occur in vacuum, to processes in many body plasma systems. When the plasma 
is in thermodynamical equilibrium one can either follow the so-called imaginary 
time formalism, see e.g. (I22|, I23 [) or the real time (Schwinger-Keldysh) formalism 
(|24[|25p . The latter can also be extended to non-equilibrium situations (|26U27p . 

The perturbative expansion expressed in terms of Feynman diagrams allows one 
a systematic computation of various quantities. However, in order to obtain a 
gauge invariant finite result, one often has to resum a class of diagrams as required 
by the Hard Loop Approach (28,29,30) (the real-time formulation is discussed in 
([3T]) ). The approach, which was first developed for equilibrium systems (28 ll29|130[) 
(for a review see ([321) ) and then extended to the non-equilibrium case (33, 3"4ll35"1) . 
distinguishes soft from hard momenta. In the case of ultrarelativitic QED plasmas 
in equilibrium, the soft momenta are of order eT while the hard momenta are of 
order T with T being the plasma temperature. One obviously assumes that l/e> 
1. The Hard Loop Approach deals with soft collective excitations generated by 
hard plasma particles which dominate the distribution functions. 

As an example, we consider the polarization tensor given by Equation El which 
was obtained within the kinetic theory in Sect. 12. ll We restrict ourselves to ultra- 
relativitic QED plasmas. In the lowest order of the perturbative expansion, the 
polarization tensor or photon self-energy is given by the diagram shown in Fig- 
ured! The tensor can be decomposed into the vacuum and medium contributions. 




Q 



(35) 
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Figure 1: The lowest order contribution to the QED polarization tensor. 

The first one requires an usual renormalization because of a ultraviolet divergence, 
whereas the medium part appears to be ultraviolet finite. One reproduces Equa- 
tion [5] applying to the diagrammatic result the Hard Loop Approximation which 
requires that the energy and momentum (a;, k) of the external photon line are 
much smaller than the momentum (p) of the electron loop. Then, it appears that 
the vacuum part can be neglected, as it is much smaller than the medium part. 
In the case of an ultrarelativistic equilibrium EMP, Equations [7] were derived 
diagrammatically in (|36|l37p . In QGP the lowest order polarization tensor (gluon 
self energy) includes one- loop diagrams with internal gluon and ghost lines. The 
final result for the gluon polarization tensor in the high-temperature approxima- 
tion essentially coincides with the QED expression. The color degrees of freedom 
enter through the trivial color factor 5 a b- In the case of equilibrium QGP, one 
additionally replaces in Equation [7J the thermal photon mass by a thermal gluon 
mass given by 

2 g 2 T 2 f,,N f \ 
m * = — v 1+ 6j' (36) 
where Nf indicates the number of light quark flavors. 

The Hard Loop Approach can be nicely formulated in terms of an effective 
action. Such an action for an equilibrium system was derived diagrammatically 
in (|29p and in the explicitly gauge invariant form in (|30l) . The equilibrium Hard 
Loop action was also found within the semiclassical kinetic theory (|38U 39D. The 
action was generalized ([33,35) for non-equilibrium systems which, however, are 
on average locally color neutral, stationary and homogeneous. 

The starting point was the effective action which describes an interaction of 
classical fields with currents induced by these fields in the plasma. The lagrangian 
density is quadratic in the gluon and quark fields and it equals 
'1 



C 2 (x) = -J &S 



2 >• 



-A*Jx)IL%(x - y)K(y) + *(*)£(* - y)0>(y) 



(37) 



where H^ b and £ is the gluon polarization tensor and the quark self-energy, 
respectively, while A a and \& denote the gluon and quark fields. Following Braaten 
and Pisarski (|30p , the lagrangian from Equation [37] was modified to comply with 
the requirement of gauge invariance. The final result, which is non-local but 
manifestly gauge invariant, is 

N 2 — 1 ~ - v ^ 
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where F£ v is the strength tensor and D denotes the covariant derivative; /(p) 
and /(p) are the effective parton distribution functions defined as /(p) = n(p) + 
n(p) + 2N c n g (p) and /(p) = n(p) + n(p) + 2n 9 (p); n(p), n(p) and n s (p) are 
the distribution functions of quarks, antiquarks and gluons of a single-color com- 
ponent in a homogeneous and stationary plasma which is locally and globally 
colourless; the spin and flavor are treated as parton internal degrees of freedom. 
The quarks and gluons are assumed to be massless. The effective action given 
by Equation [38] generates n— point functions which obey the Ward-Takahashi 
identities. Equation [38] holds under the assumption that the field amplitude is 
much smaller than T/g where T denotes the characteristic momentum of (hard) 
partons. 

3 COLLECTIVE PHENOMENA 

The most characteristic feature of the EM and QCD plasma, which results from a 
long range interaction governing both systems, is a collective behavior which leads 
to specific plasma phenomena like screening, plasma oscillations, instabilities, etc. 

Since the electromagnetic and chromodynamic polarization tensors, which are 
obtained in the linear response analysis, are essentially the same, the collective 
effects in EMP and QGP are very similar in the linear response regime. As 
our discussion is limited to this regime, mostly the electromagnetic plasma is 
considered in this section. 

3.1 Screening 

We start with screening of electric charges in the plasma. To discuss the effect, 
let us consider an electric field generated by the point-like charge q moving with a 
velocity v in the plasma. The problem is studied in numerous plasma handbooks 
e.g. in (]2j). The induction vector obeys the Maxwell equation 

V-D(x) =4vrg(5 (3) (r- vt) , 

with x = (t, r). After the Fourier transformation, which is defined by Equation [3] 
the induction vector reads 

ik ■ D(fc) = 8Tr 2 q 5(to - k • v) , (39) 

where k = (u, k). The induction vector D(/c) is related to the electric field E(fc) 
through the dielectric tensor e lJ (k) as 

D i {k) = e ij {k)Ei{k) . (40) 

We note that the dielectric tensor E l3 {k), which carries information on the elec- 
tromagnetic properties of a medium, can be expressed through the polarization 
tensor as 

= S ij + . (41) 
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In the isotropic plasma there are only two independent components of the dielec- 
tric tensor et and el which are related to e 13 as 

£ ij (k) = E T {k)(8 ij - k' l k j /k 2 ) + £ L (k) k l k j /k 2 . (42) 

Using Equations 001 021 and expressing the electric field E through the scalar 
cf) and vector A potentials (E(fc) = — i\n<j)(k) + icoA(k)) in the Coulomb gauge 
(k • A(k) = 0), one finds the electric potential in a medium (the wake potential) 

/^3£, gik(r-vt) 
(2ir) d e l (uj = v • k, k) K 2 

Let us first consider the simplest case of the potential generated by a static 
(v = 0) charge. Using Equations [H [7J £z,(0, k) of an ultrarelatvistic electron- 
positron plasma is found as 

e L (0,k) = l + ^|, (44) 

where mo is the so-called Debye mass given by m 2 D = e 2 T 2 /3 = 3m 2 = IIl(0, k). 
Then, Equation 03] gives the well-known screened potential 

0(r) = ^ e~ mDr , (45) 
r 

with r = |r|. As seen, the inverse Debye mass has the interpretation of the 
screening length of the potential. Since the average inter-particle spacing in the 
ultrarelatvistic plasma is of order T _1 , the number of particles in the Debye sphere 
(the sphere of the radius m^ 1 ) is of order e -3 which is, as already mentioned in 
the Introduction, much larger that than unity in the weakly coupled plasma 
(1/e 2 3> 1). This explains the collective behavior of the plasma, as motion of 
particles from the Debye sphere is highly correlated. 

For v the potential given by Equation 03] has a rich structure. In the 
context of QGP it has been discussed in ()4O1 I41 1021) . showing that it can exhibit 
attractive contributions even between like-sign charges in certain directions (|40p . 
For a supersonic particle, the potential can reveal a Mach cone structure associ- 
ated with Cerenkov radiation when electromagnetic properties of the plasma are 
appropriately modeled (|41U42p . 



3.2 Collective modes 

Let us consider a plasma in a homogenous, stationary state with no local charges, 
no currents. As a fluctuation or perturbation of this state, there appear local 
charges or currents generating electric and magnetic fields which in turn interact 
with charged plasma particles. Then, the plasma reveals a collective motion which 
classically is termed as plasma oscillations. Quantum-mechanically we deal we 
quasi-particle collective excitations of the plasma. 

The collective modes are found as solutions of the dispersion equation obtained 
from the equation of motion of the Fourier transformed electromagnetic potential 
A^{k) which is 

k 2 g nv _ k n k v _ nf*^ 



A v {k) = , 



(46) 
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where the polarization tensor IP 11 ' contains all dynamical information about the 
system. The general dispersion equation is then 

det \k 2 g^ - k^k u - IP" (jfe)l = . (47) 

Due to the transversality of IP* 1 {k) not all components of H^(k) are inde- 
pendent from each other, and consequently the dispersion Equation |4"TI which 
involves a determinant of a 4 x 4 matrix, can be simplified to the determinant 
of a 3 x 3 matrix. For this purpose one usually introduces the dielectric tensor 
e lJ (k) related to the polarization tensor by Equation 0TJ Then, the dispersion 
equation gets the form 

detfk 2 ^' - k { k j - uj 2 e ij {k)\ = . (48) 

The relationship between Equation [47] and Equation [48] is most easily seen in the 
Coulomb gauge when = and k • A(k) = 0. Then, E = iujA. and Equation 06] 
is immediately transformed into an equation of motion of E(fc) which further 
provides the dispersion Equation [48] 

As expressed by Equation l42l there are only two independent components of the 
dielectric tensor (er(fc) and ££,(&)) in the isotropic plasma. Then, the dispersion 
Equation [48] splits into two equations 

e T {k) = k 2 /u; 2 , e L {k) = . (49) 

Solutions of the dispersion equations w(k), with a complex, in general, fre- 
quency u, represent plasma modes which classically are, as already mentioned, the 
waves of electric and /or magnetic fields in the plasma while quantum-mechanically 
the modes are quasi-particle excitations of the plasma system. If the imaginary 
part of the mode's frequency 3u> is negative, the mode is damped - its amplitude 
exponentially decays in time as e 3aJ *. When Qlo = we have a stable mode with 
a constant amplitude. Finally, if > 0, the mode's amplitude exponentially 
grows in time - there is an instability. 

When the electric field of a mode is parallel to its wave vector k, the mode is 
called longitudinal. A mode is called transverse when the electric field is trans- 
verse to the wave vector. The Maxwell equations show that the longitudinal 
modes, also known as electric, are associated with electric charge oscillations; the 
transverse modes, also known as magnetic, are associated with electric current 
oscillations. 

The collective (boson) modes in the equilibrium ultrarelativistic plasma are 
shown in Figure [2] There are longitudinal modes also called plasmons and trans- 
verse modes. Both start at zero momentum at the plasma frequency which is 
identical to the thermal photon (or gluon) mass, lo v = m 7 . The dispersion re- 
lations lies above the light cone (u; > |k|) showing that the plasma waves are 
undamped (no Landau damping) in the high-temperature limit. As explained 
in Sect. 13.3, 1} the Landau damping, which formally arises from the imaginary 
part of the polarization tensor given by Equation [7] at lu 2 < k 2 , occurs when the 
energy of the wave is transferred to plasma particles moving with velocity equal 
to the phase velocity (w/|k|) of the wave. If the phase velocity is larger than the 
speed of light such a transfer is obviously not possible. 
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Figure 2: Dispersion relation of longitudinal and transverse plasma waves. 



3.2.1 TWO-STREAM SYSTEM As an example of the rich spectrum of 
collective modes, we consider the two-stream system within the hydrodynamic 
approach when the effect of pressure gradients is neglected. Details of the analysis 
can be found in (|2ip . The dielectric tensor provided by the polarization tensor 
from Equation [28] is 



'(*.k>-(i-£) 



4tt 



u 2 „ e a + Pa 



v l a k j + v^k 1 (uj 2 - k 2 



UJ 



k • v Q (lo — k • v c 



(50) 



where v a is the hydrodynamic velocity related to the hydrodynamic four-velocity 



u£; u>p is the plasma frequency given as 



a ~i~ Pc 



(51) 



The index a, which labels the streams and plasma components, has four values, 
a = L—,L+,R—,R+. The first character labels the stream ( l R J for right and 
'L' for left) while the second one labels the plasma component ('+' for positive 
and '— ' for negative). For simplicity we assume here that the streams are neutral 
and identical to each other and their velocities, which are chosen along the axis 
z, are opposite to each other. Then, 



n 

P 

e 



riL- 
PL- 
QL- 



PL+ = PR- 
~<1L+ = QR- 



= n R+ , 
PR+ , 
= -QR+ 



VL- 



VL+ : 



-VR- 



-VR+ 



(52) 



and the plasma frequency equals uj 2 = 167re n /(e + p). 

The wave vector is first chosen to be parallel to the axis x, k = (k,0,0). Due 
to Equations [52l the off-diagonal elements of the matrix in Equation 08] vanish 
and the dispersion equation with the dielectric tensor given by Equation [50] is 



(w 2 -Wp)(w 



k 2 )(u? 



u; - k 



X 2 



k 2 



UJ 







(53) 



18 



Mrowczynski & Thoma 



Up and a stable transverse mode with uj 2 



to 2 + k 2 . There are 



X 2 + k 2 ± 



(oj 2 - X 2 + k 2 ) 2 + 4X 2 k 2 



(54) 



where A 2 = u 2 v 2 . As solutions of the equation, one finds a stable longitudinal 
mode with to 2 
also transverse modes with 

2 1 

± 2 

As seen, uj\ > but uo 2 _ < 0. Thus, the mode uj + is stable and there are two 
modes with pure imaginary frequency corresponding to uPl < 0. The first mode 
is overdamped while the second one is the well-known unstable Weibel mode f)43|) 
leading to the filamentation instability. A physical mechanism of the instability 
is explained in Sect. 13.3.21 

The wave vector, as the stream velocities, is now chosen along the z— axis i.e. 
k = (0,0, k). Then, the matrix in Equation [48] is diagonal. With the dielectric 
tensor given by Equation [501 the dispersion equation reads 



(u, 2 



-Ul r 



-k 



2\2 



U) —uo r 



-ti- 



ki) 



; + 



(k 2 - io 2 )v 2 



kv 



+ 



(k 2 - co 2 )v 2 



kv 2(u — kv) 2 lo + kv 2(u + kv) 



= . 

(55) 

ujp + k' z . The longitudinal 



There are two transverse stable modes with uj 2 - 2 1 ' 
modes are solutions of the above equation which can be rewritten as 

1 1 



1 



(uo — kv)'' 



+ 



0, 



(56) 



where uJq = uj 2 /2^ 2 with 7 = (1 — v 2 
x = uj/ojq, y = kv/ujQ, Equation [56] is 



(u + kv) 2 _ 
l l 2 . With the dimensionless quantities 



(x 2 - y 2 ) 2 - 2x 2 -2y 2 = 0, 



(57) 



and it is solved by 



r 2 

X± 



y~ + 1 ± v /4y 2 + 1 . (58) 

As seen, x\ is always positive and thus, it gives two real (stable) modes; x 2 _ is 
negative for < y < \[2 and then, there are two pure imaginary modes. The 
unstable one corresponds to the two-stream electrostatic instability. A physical 
mechanism of the instability is explained in Sect. 13.3.11 



3.3 Instabilities 

Presence of unstable modes in a plasma system crucially influences its dynamics. 
Huge difficulties encountered by the half-a-century program to built a thermonu- 
clear reactor are just related to various instabilities experienced by a plasma 
which make the system's behavior very turbulent, hard to predict and hard to 
control. 

There is a large variety of instabilities - the history of plasma physics is said 
to be a history of discoveries of new and new instabilities. Plasma instabilities 
can be divided into two general groups (1) hydrodynamic instabilities, caused by 
coordinate space inhomogeneities, (2) kinetic instabilities due to non-equilibrium 
momentum distribution of plasma particles. 
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Figure 3: The mechanism of energy transfer between particles and fields. 

The hydrodynamic instabilities are usually associated with phenomena occur- 
ring at the plasma boundaries. In the case of QGP, this is the domain of highly 
non-perturbative QCD where nonAbelian nature of the theory is of crucial im- 
portance. Then, the behavior of QGP is presumably very different from that of 
EMP, and thus, we will not speculate about possible analogies. 

The kinetic instabilities are simply the collective modes with positive Qco intro- 
duced in Sect. 13.21 and found in Sect. 13.2.11 in the specific case of the two-stream 
system. Thus, we have longitudinal (electric) and transverse (magnetic) insta- 
bilities. In the non-relativistic plasma the electric instabilities are usually much 
more important than the magnetic ones, as the magnetic effects are suppressed by 
the factor v 2 /c 2 where v is the particle's velocity. In the relativistic plasma both 
types of instabilities are of similar strength. As will be discussed later on, the 
electric instabilities occur when the momentum distribution of plasma particles 
has more than one maximum, as in the two-stream system. A sufficient condi- 
tion for the magnetic instabilities appears to be an anisotropy of the momentum 
distribution. 

3.3.1 MECHANISM OF ELECTRIC INSTABILITY Let us consider 
a plane wave of the electric field with the wave vector along the z axis. For 
a charged particle, which moves with a velocity v = p z /E p equal to the phase 
velocity of the wave = u/k, the electric field does not oscillates but it is 
constant. The particle is then either accelerated or decelerated depending on the 
field's phase. For an electron with v = chances to be accelerated and to be 
decelerated are equal to each other, as the time intervals spent by the particle in 
the acceleration zone and in the deceleration zone are equal to each other. 

Let us now consider electrons with the velocities somewhat smaller the phase 
velocity of the wave. Such particles spend more time in the acceleration zone 
than in the deceleration zone, and the net result is that the particles with v < 
are accelerated. Consequently, the energy is transferred from the electric field 
to the particles. The particles with v > v a spend more time in the deceleration 
zone than in the acceleration zone, and thus they are effectively decelerated - the 
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Figure 4: Momentum distribution of the plasma-beam system. 



energy is transferred from the particles to the field. If the momentum distribution 
is such that there are more electrons in a system with v < v $ than with v > tu, 
the wave looses energy which is gained by the particles, as shown in the left- 
hand-side of Figure This is the mechanism of famous collisionless Landau 
damping of the plasma oscillations. If there are more particles with v > than 
with v < Vff,, the particles loose energy which is gained by the wave, as in the 
right-hand-side of Figure [3j Consequently, the wave amplitude grows. This is the 
mechanism of electric instability. As explained above, it requires the existence 
of the momentum interval where f n (p) grows with p. Such an interval appears 
when the momentum distribution has more than one maximum. This happens 
in the two-stream system discussed in Sect. 13.2.11 or in the system of a plasma 
and a beam shown in Figured! 

3.3.2 MECHANISM OF MAGNETIC INSTABILITY Since the mag- 
netic instabilities appear to be relevant for QGP produced in relativistic heavy-ion 
collisions (see below), we discuss them in more detail. Let us first explain follow- 
ing (|44p . how the unstable transverse modes are initiated. For this purpose we 
consider a plasma system which is homogeneous but the momentum distribution 
of particles is not of the equilibrium form, it is not isotropic. The system is on av- 
erage locally neutral (0' M (x)) = 0) but current fluctuations are possible, and thus 
the correlator {j^{x\)j 1 ' (x 2 )) is i n general non-zero. Since the plasma is assumed 
to be weakly coupled, the correlator can be estimated neglecting the interaction 
entirely. Then, when the effects of quantum statistics are also neglected, the 
correlator equals 

M^(t,x) d ^ f (^(t 1 ,x 1 )/(t 2) x 2 ))=^^ J _0_ ^ /n(p) 5(3) (x _ v i ))(59) 

where v = p/E p and (t,x) = (i 2 — ii,x 2 — xi). Due to the average space-time 
homogeneity, the correlator given by Equation [59] depends only on the difference 
(i 2 — ti, x 2 — xi). The space-time points (ii, xi) and (t 2 , x 2 ) are correlated in the 
system of non-interacting particles if a particle travels from (ti,xi) to (t 2 ,x 2 ). 
For this reason the delta 8^ (x — vt) is present in Equation [59l The sum and 
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Figure 5: The mechanism of filamentation instability, see text for a description. 

momentum integral represent summation over all particles in the system. The 
fluctuation spectrum is found as the Fourier transform of Equation 1591 i.e. 

M k) = £ ql J ^0 /n ( p ) 2n5(u - kv) . (60) 

To further study the fluctuation spectrum, the particle's momentum distri- 
bution has to be specified. We will present here only a qualitative discussion of 
Equations 1591 IBTil assuming that the momentum distribution is strongly elongated 
in one direction which is chosen to be along the z axis. Then, the correlator M zz 
is larger than M xx or M yy . It is also clear that M zz is the largest when the wave 
vector k is along the direction of the momentum deficit, as in such a case the delta 
function 5(oj — kv) does not much constrain the integral in Equation 1601 Since the 
momentum distribution is elongated in the z direction, the current fluctuations 
are the largest when the wave vector k is the x — y plane. Thus, we conclude 
that some fluctuations in the anisotropic system are large, much larger than in 
the isotropic one. An anisotropic system has a natural tendency to split into 
the current filaments parallel to the direction of the momentum surplus. These 
currents are seeds of the transverse unstable mode known as the filamentation or 
Weibel instability (133]) which was found in the two-stream system discussed in 
Sect. I5XT1 

Let us now explain in terms of elementary physics why the fluctuating currents, 
which flow in the direction of the momentum surplus, can grow in time. The form 
of the fluctuating current is chosen to be 

j(x) = j e z cos(k x x) , (61) 

where e z is the unit vector in the z direction. As seen in Equation [BTl there are 
current filaments of the thickness vr/l^l with the current flowing in the opposite 
directions in the neighboring filaments. The magnetic field generated by the 
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current from Equation [61] is given as 

j 

B(x) = An — e y s'm(k x x) , 

and the Lorentz force acting on the particles, which fly along the z direction, 
equals 

j 

F(x) = q v x B(x) = —Airq v z — e x sm(k x x) , 

kx 

where q is the particle's electric charge. One observes, see Figure that the force 
distributes the particles in such a way that those, which positively contribute to 
the current in a given filament, are focused in the filament center while those, 
which negatively contribute, are moved to the neighboring one. Thus, the initial 
current is growing and the magnetic field generated by this current is growing as 
well. The instability is driven by the the energy transferred from the particles 
to fields. More specifically, the kinetic energy related to the motion along the 
direction of the momentum surplus is used to generate the magnetic field. 
3.3.3 ROLE OF INSTABILITIES As already mentioned, there is a large 
variety of plasma instabilities which strongly influence numerous plasma char- 
acteristics. Not much is known of the hydrodynamic instabilities of QGP and 
if they exist, they belong to the highly non-perturbative sector of QCD which 
is still poorly understood. As explained in Sect. 13.3. H the electric instabilities 
occur in a two-stream system, or more generally, in systems with a momentum 
distribution of more than one maximum. While such a distribution is common 
in EMP, it is rather irrelevant for QGP produced in relativistic heavy- ion col- 
lisions where the global as well as local momentum distribution is expected to 
monotonously decrease in every direction from the maximum. The electric in- 
stabilities are absent in such a system but a magnetic unstable mode, which has 
been discussed in Sect. 13.3.21 is possible. The filamentaion instability was first 
argued to be relevant for the QGP produced in relativistic heavy-ion collisions 
in (|45U46^[44|) . A characteristic time of instability growth was estimated (|45p46p 
to be shorter or at least comparable to other time scales of the parton system 
evolution. The mechanism of instability growth was also clarified (|44p . The early 
arguments were substantiated in the forthcoming analytic calculations (47, 4"8jl49p 
and numerical simulations (50,51,52,53,54). 

A main consequence of instabilities is a fast equilibration of the weakly coupled 
plasma. The problem is of particular interest as the experimental data on heavy- 
ion collisions, where QGP production is expected, suggest that an equilibration 
time of the parton system is below 1 fm/c (|55l) , A whole scenario of instabilities 
driven equilibration is reviewed in ()56j) . Here we only mention the main points, 
starting with an observation that collisions of charged particles are not very 
efficient in redistributing particle momenta, as the Rutherford cross section is 
strongly peaked at a small momentum transfer. One either needs many frequent 
but soft collisions or a few rare but hard collisions to substantially change a 
particle's momentum. As a result the inverse time of collisional equilibration 
of QGP is of order g 4 ln(l/g)T (|16j) where T is the characteristic momentum of 
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quarks or gluons. It appears that the momentum distribution is isotropized due 
to instabilities within the inverse time of order gT (|57p . If 1/g 3> 1, the collisional 
equilibration is obviously much slower. As discussed in Chapter U the situation 
changes in strongly-coupled plasmas. 

When the instabilities grow, the system becomes more and more isotropic 
(|46U57p because the Lorentz force changes the particle's momenta and the growing 
fields carry an extra momentum. To explain the mechanism let us assume that 
initially there is a momentum surplus in the z direction. The fluctuating current 
flows in the z direction with the wave vector pointing in the x direction. Since the 
magnetic field has a y component, the Lorentz force, which acts on partons flying 
along the z axis, pushes the partons in the x direction where there is a momentum 
deficit. The numerical simulation ([52]) shows that growth of the instabilities is 
indeed accompanied with the system's fast isotropization. 

The system isotropizes not only due to the effect of the Lorentz force but also 
due to the momentum carried by the growing field. When the magnetic and 
electric fields are oriented along the y and z axes, respectively, the Poynting 
vector points in the direction x that is along the wave vector. Thus, the momen- 
tum carried by the fields is oriented in the direction of the momentum deficit of 
particles. 

Although the scenario of instabilities driven equilibration looks very promising, 
the problem of thermalization of QGP produced in heavy-ion collisions is far from 
being settled. It has been shown (|58p that inter-parton collisions, which have been 
modeled using the BGK collision term ([10]), reduce the growth of instabilities and 
thus slow down the process of equilibration. The equilibration is also slowed down 
due to expansion of QGP into vacuum (|59ll60j) which is a characteristic feature 
of QGP produced in relativistic heavy-ion collisions. Finally, the late stage of 
instability development, when nonAbelian effects are crucially important, appears 
to be very complex (61,62) and it is far from being understood. 

As already mentioned, instabilities influence various plasma characteristics. 
In particular, it is known (|63p that turbulent magnetic fields generated in the 
systems, which are unstable with respect to transverse modes, are responsible for 
a reduction of plasma viscosity. Then, an anomalously small viscosity, which is 
usually associated with strongly coupled systems, can occur in weakly coupled 
plasmas as well. Recently, it has been argued (I64^l65p that the mechanism of 
viscosity reduction is operative in the unstable QGP. 

3.4 Energy loss 

A charged particle which moves across the plasma changes its energy due to sev- 
eral processes (2). When the particle's energy (E) is comparable to the plasma 
temperature (T), the particle can gain energy due to interactions with field fluc- 
tuations. (In context of QGP the problem was studied in ([68).) A fast particle 
with E ^> T looses its energy and main contributions come from collisions with 
other plasma particles and from radiation. In the following we discuss the energy 
loss of a fast particle, as the problem is closely related to jet quenching which was 
suggested long ago as a signature for the QGP formation in relativistic heavy-ion 
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collisions ([66^ 167 

Let us start with the collisional energy loss. The particle's collisions are split 
into two classes: hard with high-momentum transfer, corresponding to the col- 
lisions with plasma particles and, soft with low-momentum transfer dominated 
by the interactions with plasma collective modes. The momentum is called 'soft' 
when it is of order of the Debye mass, mo, or smaller and it is 'hard' when it is 
larger than mrj. 

The soft contribution to the energy loss, which can be treated in a classical 
way, is often called the process of plasma polarization. It leads to the energy loss 
per unit time given by the formula 




soft 



Jd 3 xj(x)E(x) , (62) 



where E is the electric field induced in the plasma by the particle's current j 
which is of the form = qv8^(x — vt). The field can be calculated by means 
of the Maxwell equations. After eliminating the magnetic field, one finds the 
equation 




k*hP 



Ej (k) = -f(k) 



Since we consider equilibrium plasma, which is isotropic, one introduces the lon- 
gitudinal (el) and transverse (el) components of e 13 . Then, Equation E2] is 
manipulated to 

fdE\ 47rie 2 f d 3 k f u v 2 - co 2 /k 2 \ 

\dx~) soit ~ ~J j2^{k 2 e L (k) + uj[e T {k) - k 2 /u; 2 ] J ' ^ 

which gives the energy loss per unit length. This formula describes the effect of 
medium polarization. However, three comments are in order here. 

1. Equation [63] includes the charge self- interaction signaled by the ultraviolet 
divergence of the integral from Equation 1631 The self-interaction is removed 
by subtracting from Equation [63] the vacuum expression with el = ex = 1. 

2. Poles of the function under the integral Equation [63] correspond to the 
plasma collective modes as given by the dispersion Equations 09] Therefore, 
the explicit expressions of el and el are not actually needed to compute the 
integral in Equation [63] The knowledge of the spectrum of quasiparticles 
appears to be sufficient. 

3. Equation [63] is derived in the classical approximation which breaks down for 
a sufficiently large k. Therefore, an upper cut-off is needed. The interaction 
with k above the cut-off, which, as already mentioned, is of order of the 
Debye mass, should be treated as hard collisions with plasma particles. 



The energy loss per unit length due to hard collisions is 

? J ni{k) [flux factor] J dn ^n u ' (64) 




hard 
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where the sum runs over particle species distributed according to rii(k), v = E—E' 
is the energy transfer, and do % jd£l is the respective differential cross section. 

Combining Equations 1631 and 1641 one finds the complete collisional energy loss. 
The calculations of the energy loss of a fast parton in the QGP along the lines 
presented above were performed in ([69t l70 |) . Systematic calculations of the col- 
lisional energy loss using the Hard Thermal Loop resummation technique were 
given in (|72U7ip with a result which is infrared finite, gauge invariant, and com- 
plete to leading order. Recently, the calculations of the collisional energy loss 
have been extended to anisotropic QGP (73). 

It was realized that a sizeable contribution to the quark's energy loss comes 
from radiative processes (|74p . The problem, however, appeared to be very com- 
plex because quark's successive interactions in the plasma cannot be treated as 
independent from each other and there is a destructive interference of radiated 
gluons known as the Landau-Pomeranchuk-Migdal effect (|75p . There are numer- 
ous papers devoted to the radiative energy loss and the whole problem is reviewed 
in ([76]) . A general conclusion of these studies is that the energy lost by a fast 
light quark quadraticaly (not linearly) depends on the path traversed in QGP, as 
the radiative energy loss dominates over the collisional. Recent experiments at 
RHIC show (|77ll78p . however, that heavy quarks, whose radiative energy loss is 
significantly suppressed, are strongly dragged in the QGP medium. It may sug- 
gest that the collisional energy loss should be actually enhanced as theoretically 
predicted in ([79]) . 

4 STRONGLY COUPLED PLASMAS 

Our discussion of the collective phenomena presented in Sect. [3] was limited to the 
weakly interacting plasmas with the coupling constant much smaller than unity. 
However, QGP produced in ultrarelativistic heavy-ion collisions is presumably 
strongly coupled (sQGP) as the temperature is never much larger than Aqcd> 
and the regime of asymptotic freedom is not reached. QGP is certainly a strongly 
interacting system close to the confinement phase transition. There are indeed 
hints in the extensive experimental material collected at RHIC (80,81,82,83) that 
the matter produced at the early stage of nucleus-nucleus collisions RHIC is in 
the form of sQGP for a few fm/c. In particular, the characteristics of elliptic flow 
and particle spectra, which are described well by ideal hydrodynamics, seem to 
indicate a fast thermalization and small viscosity of the plasma. Both features are 
naturally explained assuming a strong coupling of the plasma (55,84,85,86,87). 

Although a fast thermalization (|56p as well as a small viscosity (|64i|65l) can also 
be explained by instabilities, the idea of sQGP certainly needs to be examined. 
However, the theoretical tools presented in Chapter [2] implicitly or explicitly 
assume small coupling constant and they are of limited applicability. A powerful 
approach, which can be used to study sQGP is the lattice formulation of QCD, for 
a review see (|88p . However, lattice QCD calculations, which are mostly numerical, 
encounter serious problems to incorporate quark degrees of freedom. It is also 
very difficult to analyze time dependent plasma characteristics. 
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Strongly coupled conformal field theories such as supersymmetric QCD can be 
studied by means of the so-called AdS/CFT duality ([89]) . Although some very 
interesting results on the conformal QGP were obtained in this way, see e.g. (|90p 
and references therein, relevance of these results for QGP governed by QCD, not 
by supersymmetric QCD, is unclear. Thus, the question arises what we can learn 
about sQGP from strongly coupled EMP. 

We first note that most of EMP in nature and technological applications are 
weakly coupled i.e. the interaction energy between the plasma particles is much 
smaller than their thermal (kinetic) energy. This is because strongly coupled 
plasmas require a high particle density and/or low temperature, at which usu- 
ally strong recombination occurs and the plasma state vanishes. Exceptions are 
the ion component in white dwarfs, metallic hydrogen and other states of dense 
warm matter in the interior of giant planets, short-living dense plasmas produced 
by intense laser or heavy ion beams or in explosive shock tubes, dusty (or com- 
plex) plasmas, and two-dimensional electron systems on liquid helium (19111920930 . 
Therefore, it is a real challenge to study strongly coupled EMP both theoretically 
and experimentally. 

In nonrelativistic EMP the interaction energy is given by the (screened) Coulomb 
potential. The Coulomb coupling parameter defined by 

Q 2 

r = ±= (65) 

a 1 

distinguishes between weakly coupled, r « 1, and strongly coupled plasmas, 
r~l. Here q is the particle charge, a the interparticle distance and T the kinetic 
temperature of the plasma component (electrons, ions, charged dust grains) under 
consideration. In the case of a degenerate plasma, e.g. the electron component 
in a white dwarf, the kinetic energy T is replaced by the Fermi energy. Due to 
the strong interaction, the plasma can behave either as a gas or a liquid or even 
a solid (crystalline) system. 

The case of an one-component plasma (OCP) with a pure Coulomb interaction 
(a single species of charged particles in an uniform, neutralizing background) has 
been studied as a reference model for strongly coupled plasmas using simple mod- 
els as well as numerical simulations in great detail (|91 j) . For T > 172 the plasma 
was shown to form regular structures (Coulomb crystallization) (I94p . Below this 
critical value the OCP is in the supercritical state. For values of T larger than 
about 50, it behaves like an ordinary liquid, while for small values below unity 
like a gas. Only if T is large enough, the usual liquid behavior (Arrhenius law 
for the viscosity, Stokes-Einstein relation between self-diffusion and shear viscos- 
ity, etc.) appears due to caging of the particles (a single particle is trapped for 
some period of time in the cage formed by its nearest neighbors). For values 
of r, which are smaller than about 50, caging is not sufficiently strong and the 
system shows complicated, not yet understood transport properties. However, 
the short-range ordering typical for liquids shows up already for T > 3 (95). A 
gas-liquid transition, requiring a long-range attraction and a short-range repul- 
sion, e.g. Lennard-Jones potential, does not exist in the OCP with particles of 
like-sign charges. 
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Figure 6: Phase diagram of a strongly coupled Yukawa system from Ref. ([96 



In realistic systems with a screened Coulomb interaction (Yukawa potential), 
the phase diagram can be shown in the T-K-plane, where k = cl/Xd is the distance 
parameter with Xr> being the Debye screening length. Numerical simulations 
based on molecular dynamics lead to the phase diagram shown in Figure [6] (96). 

The first quantity of interest of sQGP is the coupling parameter. In analogy 
to nonrelativistic EMP, it is defined as ([9j 



1.5 - 5, 



(66) 



4vraT 

where C is the Casimir invariant (C = 4/3 for quarks and C = 3 for gluons), 
a ~ 0.5 fm is the interparticle distance, and T ~ 200 MeV is the QGP temper- 
ature corresponding to a strong coupling constant g ~ 2. The factor 2 in the 
numerator comes from taking into account the magnetic interaction in addition 
to the static electric (Coulomb) interaction, which are of the same magnitude 
in ultrarelativistic plasmas. The factor 4-7T in the denominator comes from using 
the Heavyside-Lorentz system in QCD as discussed in the Introduction. The 
distance parameter k of the QGP under the above conditions is rather small, 
typically between 1 and 3 ([98]) . 

It should be noted, that we have assumed here a classical interaction potential 
corresponding to a one-gluon exchange. However, an effective potential taking 
into account higher order and non-perturbative effects may be much larger. This 
might be related to the fact that experimental data suggest a cross section en- 
hancement for the parton interaction by more than an order of magnitude (see 
below). Hence the effective coupling parameter might be up to an order of mag- 
nitude larger than Equation 1661 

As discussed above, comparison to the OCP model as well as experimental 
data suggest that QGP close to the confinement phase transition could be in a 
liquid phase. So, the question arises whether there is a phase transition from 
a liquid to a gaseous QGP, as sketched in Figure For such a transition, a 
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Figure 7: Sketch of a phase diagram of strongly interacting matter with a possible 
gas-liquid transition in the QGP phase. T denotes the temperature while p is the 
baryon density. 

Lennard-Jones type interaction between the partons is required. However, the 
parton interaction in perturbative QCD is either purely repulsive or attractive in 
the various interaction channels, e.g quark-antiquark or diquark channel. Due to 
non-linear effects caused by the strong coupling, however, attractive interactions 
can arise even in the case like-sign charges (see e.g. (99)) leading to Lennard- 
Jones type potentials. Hence a gas-liquid transition in QGP with a critical point, 
proposed in (|87p . is not excluded and deserves further investigation. 

An important quantity, which is very useful in theoretical and experimental 
studies of strongly coupled systems on the microscopic level, in particular in 
fluid physics (100), are correlation functions. In particular, the pair correlation 
function and the static structure function provide valuable information on the 
equation of state of the system (jlOOp . The extension of the approach to the QGP 
has been proposed in ([98]) . 

The static density-density autocorrelation function is defined for a classical 
system as (fT00ll9Tj) 

G(r) = ^ Jd 3 r' (p(r + v',t)p(r',t)), 
where N is the total number of particles and 

N 
i=i 

is the local density of point particles with Ti(t) denoting the position of i— th 
particle at time t. The density-density autocorrelation function is related to the 
pair correlation function, which is defined as 

1 N 

9(r) = x( E ^(r + r.-r,)), 
by the relation G(r) = g(r) + 5^ 3 ^(r). The static structure function, defined by 

S( P ) = ^ (p(p)p(-p)) 




Figure 8: Sketch of the static structure functions vs. momentum in the gas and 
liquid phase in arbitrary units. 



with the Fourier transformed particle density, 

p(p) = J tfVp(r) e"** 
is the Fourier transform of the density-density autocorrelation function 



S(p) = J d 3 r e- ip r G{r) 



The static structure function S(p) is constant for p / for uncorrelated par- 
ticles ([2]). The typical behavior of the function in an interacting gas and in a 
liquid is sketched in Figure [8] The oscillatory behavior is caused by short range 
correlations, corresponding to a short range ordering typical for liquids. In the 
case of OCP, the oscillations appear for T > 3 indicating a liquid behavior (albeit 
with non-standard transport properties) of the supercritical phase already for 
rather low values of T (|95|) . 

The static quark structure function can be related to the longitudinal gluon 
polarization tensor containing only the quark loop (Figured]) via 

S(p) = ~2~ f duj l (uj, p) coth-^-, 

irg z n Jo 21 

where n = N/V = (p(r)) is the average particle density in a homogeneous system. 
As a reference for the strong coupling regime, the static quark structure function 
has been calculated in the weak coupling limit by resumming the polarization 
tensor in the high-temperature limit (Equation EJ) leading to (|98j) 



2N f T 3 p 2 

where m 2 D = Nfg 2 T 2 /6 is the quark contribution to the Debye screening mass. 
The static structure function given by Equation [67] starts at zero for |p| = and 
saturates at the uncorrelated structure function S(p) = 2NfT 3 /n for large |p|. 
Such a structure function corresponds to an interacting Yukawa system in the 
gas phase (see Figure [8]). Indeed, the pair correlation function, following from 
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the Fourier transform of S(p) — 1, is 

. N f T 3 m 2 D _ mor 

2im r 

and it reproduces the Yukawa potential. 

In order to compute the structure function in strongly-coupled plasmas, molec- 
ular dynamics is used (91). Although QGP is not a classical system, as the ther- 
mal de Broglie wave length is of the same order as the interparticle distance, 
molecular dynamics might be useful as a first estimate (|87j) . Using molecular 
dynamics for a classical sQGP (|1011ll02|ll03p . the expected behavior described 
above has been qualitatively verified (llOip . In strongly coupled dense matter, 
where quantum effects are important, quantum molecular dynamics based on a 
combination of classical molecular dynamics with the density functional theory 
has been successfully applied (|104J) . A generalization to the relativistic QGP has 
not been attempted so far. As an ultimate choice, lattice QCD could be used to 
calculate the structure or pair correlation functions. This would provide a test 
for the state of sQGP as well as the importance of quantum effects by comparing 
lattice results to classical molecular simulations. 

As a last application, we consider the influence of strong coupling on the 
cross sections entering transport coefficients (shear viscosity), stopping power, 
and other dynamical quantities of the plasma. Beside higher-order and non- 
perturbative quantum effects, there is already a cross section enhancement on 
the classical level. The reason is that the Coulomb radius, defined as r c = q 2 /E 
with the particle energy E, is of the order of the Debye screening length or r c is 
even larger than in a strongly coupled plasma. Hence the standard Coulomb 
scattering formula has to be modified since the interaction with particles out- 
side of the Debye sphere contributes significantly, and consequently the inverse 
screening length cannot be used as an infrared cutoff. This modification leads, for 
example, to the experimentally observed enhancement of the so-called ion drag 
force in complex plasmas which is caused by the ion-dust interaction (105). 

In the QGP at T ~ 200 MeV, the ratio r c /A£> equals 1-5. It might enhance 
a parton cross section by a factor of 2 - 9 (97) compared to perturbative results. 
An enhanced cross section reduces the mean free path A, and consequently it 
reduces the viscosity r\ as r/ ~ A. An enhancement of the elastic parton cross 
section by more than an order of magnitude compared to perturbative results 
also explains the elliptic flow and particle spectra observed at RHIC (106|. An 
infrared cutoff smaller than the Debye mass gives a natural explanation for this 
enhancement. We also note that if the cross section is enhanced, the collisional 
energy loss grows. On the other hand, the radiative energy loss is expected to 
be suppressed in the strongly coupled QGP by the Landau-Pomeranchuk-Migdal 
effect £75]). 

Finally, we mention two examples of strongly coupled systems which have not 
been considered in QGP physics but may be of relevance. Strongly coupled 
plasmas such as two-dimensional Yukawa liquids (|107j) and dusty plasmas are 
non-Newtonian fluids, i.e. the shear viscosity depends on the shear rate (flow 
velocity) as it is well known in daily life from ketchup (shear thinning). The 
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second example concerns nanofluidics. The expanding fireball in ultrarelativistic 
heavy-ion collisions has a transverse dimension of about 20 inter-particle distances 
(about 10 fm). Fluids consisting of such a low number of layers exhibit properties 
different from large fluid systems. For example, the shear flow does not show an 
continuous velocity gradient but jumps due to the adhesive forces between two 
layers. Such an behavior has been observed for example in complex plasmas. 
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